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C3 , I. INTRODUCTION 

The observed acceleration of cosmic fluid [l| has, up to now, no definitive explanation despite of the fact that this 
discovery has dramatic implications not only for cosmology, but also for fundamental physics. According to the results 
of WMAP, PLANCK and other experiments mapping anisotropics of the cosmic microwave background, if General 
Relativity (GR) gives the correct description of gravitational interaction at cosmological level, then approximately 
O ■ 75% of cosmic energy content is an exotic, invisible, and unclustered form of dark energy. Furthermore, in order to 
explain clustering of large scale structure and stability of galactic systems, we need almost 25% of another mysterious 
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^ ' component called dark matter. No final evidence for both these dark ingredients has been found at quantum level 



O 



X 



and this status of art brought to the formulation of the so-called ACDM model which addressing the coarse-grained 
problem of cosmology does not satisfy the issue to explain the fundamental nature of dark energy and dark matter 



^ ' An alternative approach consists in dispensing with the mysterious dark energy and dark matter and asking for the 
] fact that gravitational interaction could act in different way at large scales. Historically, this approach was pursued in 
i: explaining phenomena as the precession of Mercury's perihelion which is not due to an unseen mass but to Einstein's 
] modification of Newtonian gravity. The paradigm consists in not imposing a priori the form of gravitational action 
• , but assuming that it could be dictated by data, observations and experiments at various scales. The simplest example 
from this point of view is the so-called /(_R)-gravity where "modified" gravity consists of modifications of GR that could 
become important at low curvatures, late in the matter era (infrared limit) or at high curvatures toward inflationary 
epoch (ultraviolet limit). Specifically, the Einstein-Hilbert action is linear in the Ricci scalar R and could work only 
at local (e.g. Solar System) scales. Here f{R) is a generic, non-linear function of i? 

In principle, an effective interaction action can contain several degrees of freedom: tensor, vector, and scalar, 
massless or massive. In GR only the massless spin 2 graviton propagates. When the Einstein-Hilbert action is 
^ , modified, other degrees of freedom appear. The change R — > f{R) in the action brings to a further massive scalar 
mode which can drive both the today observed cosmic acceleration and the primordial infiation, depending on its 
mass. This is analogous to the inflaton field driving the accelerated expansion of the early universe, although at a 
much lower energy scale. 

If terms quadratic in the Ricci and Riemann tensor, and possibly other curvature invariants, are included in the 
gravitational Lagrangian, / (i?, RabR"'^, Rabcd.R'^^'^'^ , ■■■ ) , massive gravitons and vector degrees of freedom appear 
In the following we restrict ourselves to f{R) theories of gravity and focus on their initial value problem ^ . For a more 
comprehensive discussion on Extended Theories of Gravity, we refer the reader to 0, 3 • 

We adopt a conservative point of view and regard /(i?)-gravity more as a toy model than the final theory of gravity, 
i.e., we consider these theories as a proof of principle that modifying gravity is a viable alternative to dark energy 
and dark matter. However, we do not feel that one can claim that any of the f{R) models proposed thus far is the 
"correct" one, or has exceptional support from the observational data. While it is true that many /(i?) models pass 
all the available experimental tests and fit the cosmological data, the same is true for many dark energy models, and 



^ f{R) gravity, with an explicit coupling to matter to R 0|, can be also used as a valid alternative to dark matter P, 
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it is currently impossible to use today available observational data to definitely discriminate between most of them, 
and between dark energy and modified gravity models. 

Modifying gravity is risky: unwanted consequences may be the violations of experimental limits on the parametrized- 
post-Newtonian (PPN) parameters at terrestrial and Solar System scales jioi, instabilities, ghosts and, as in any newly 
proposed theory, the Cauchy problem is a very delicate issue. 

f{R) gravity has a long history: its origins can be traced back to Weyl's theory in which a term quadratic in the 
Weyl tensor was added to the Einstein-Hilbert Lagrangian [TTI|. Later, f{R) gravity received the attention by many 
authors, including Eddington, Bach, Lanczos, Schrodinger, and Buchdahl. In the 1960's and 1970's, it was found that 
quadratic corrections to the Hilbert-Einstein action were necessary to improve the renormalizability of GR and 
in 1980 quadratic corrections were found to fuel inflation without the need for scalar fields |13j. Non-linear corrections 
are also motivated by string theories We refer the reader to [l^ for a comprehensive historical review. 

The prototype of f{R) gravity is the model 

f{R) = R + aR"^ - fi'^/R, (I.l) 

where a is a mass parameter related to the so called scalaron [l^, [l^ , and is a mass scale of the order of the present 
value of the Hubble parameter /i ~ _ff o ~ 10"'^^ eV. Although ruled out by its weak-field limit fid] and by a violent 
instability [17], this model gives the idea underlying modified gravity: the 1/R correction is negligible in comparison 
with R^ at the high curvatures of the early universe, and kicks in only as i? — > 0, late in the history of the universe. 
GR is restored when the i?-term is leading. 

Several models have been proposed in literature: here we discuss only general features of /(_R)-gravity addressing, 
in particular, the Cauchy problem for metric and metric-affine formulations of these theories. 

To be more specific, any theory of physics is "physically" viable if an appropriate initial value problem is suitably 
formulated. This means that, starting from the assignment of suitable initial data, the subsequent dynamical evolution 
of the physical system is uniquely determined. In this case, the problem is said well-formulated. However, also if the 
initial value problem is well-formulated, we need other properties that a viable theory has to satisfy. First of all, small 
changes and perturbations in the initial data have to produce small perturbations in the subsequent dynamics over all 
the space-time where it is defined. This means that the theory should be "stable" in order to be "predictive". Besides, 
changes in the initial data region have to preserve the causal structure of the theory. If both these requirements are 
satisfied, the initial value problem of the theory is also well-posed. 

GR has a well-formulated and well-posed initial value problem but, as for other relativistic field theories, we need 
initial value constraints and gauge choices in order to make Einstein's field equations suitable for a correct formulation 
of the Cauchy problem. The consequence of this well-posedness is that GR is a "stable" theory with a robust causal 
structure where singularities can be classified [l^ [3 • 

In this paper, we focus attention on the well-posedness of the initial value problem of metric and metric-affine 
/(i?)-theories of gravity. The aim is to prove that the Cauchy problem is well posed in different and physically 
important cases. 

The layout of the paper is the following. In Sec. II, we give the generalities of metric formulation of /(i?)-gravity 
stressing, in particular, how it can be recast in terms of scalar-tensor gravity (the so called O'Hanlon theory). After 
we discuss the Cauchy problem showing that it is both well formulated and well posed. The metric-affine formulation 
is discussed in Sec. HI. In this case, the Cauchy problem has some subtleties that have to be faced with care. While 
the well formulation is in general always valid (at least for standard matter sources [l^), the well-posedness may 
depends on the specific form of the function f{R). The relevant cases of perfect fluid and scalar field as sources are 
considered in details. Sec. IV is devoted to the discussion of the results and the conclusions. 



II. METRIC /(i?)-GRAVITY 



A. Generalities 



The action of f{R)- gravity in metric formalism is given by 

Ag) = I ^\[f{R)+Cra] ds (II. 1) 



where f{R) is a real function of the scalar curvature R{g) — g^-' Rij (with Rij :— R'\fij) and Cmds is the material 
Lagrangian. The field equations are obtained by varying pi.ip with respect to the metric gij. This yields 

nR)R^J - lfiR)9^J - S,, + V.V,/'(i?) - g^jgP'^VpVJ'iR) (II.2) 
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where E^^ 
like form 



1 (5/:„ 



is the energy-momentum tensor of perfect-fluid matter. Eqs. (jll.2[) can be recast in Einstein- 



{f{R)R^f{R))g,, 



(11.3) 



/(i?)-gravity can be re-interpreted as an O'Hanlon scalar-tensor theory, by introducing a suitable scalar field if which 
non-minimally couples with the gravity sector To see this point, let us take into account an O'Hanlon Lagrangian 

of the form 

(11.4) 



A{g, v) = j VWW^R- V{^) + C,rr] ds 



where V{ip) is the self-interaction potential. Field equations are derived by varying pi.4p with respect to both gij 
and Lp which now represents a new dynamical variable. One obtains 



Rij - 2^91] - — 



1 



(IL5a) 



R-'^^O 
dip 

We notice that, replacing pi.5bp into the trace of (|II.5ap . we obtain the scalar equation 



dip 



(II.5b) 



(11.6) 



System (jll.5|) is then equivalent to eqs. pi.5a[) together with pi.6p . Given the function f{R) in pi.l|) . we shall suppose 
that its first derivative f'{R) is invertible. In such a circumstance, it is easily seen that metric /(i?)-theories of gravity 
can be mapped onto O'Hanlon theories and vice- versa. Indeed, defined the scalar field 



and the potential 



^ := f'iR) 



Vip) :=/' [R{ip))Ri^)^f{R{ip)] 



(II.7a) 



(II.7b) 



It is straightforward to verify that, under the above hypothesis f"{R) ^ 0, eq. pi.7a[) expresses the inverse relation 
of pl.Sbp . namely 



dip 



V f'iR) 



(H., 



being the potential V{ip) defined by pi.7b[) . A direct comparison of eqs. pi.3|) with eqs. pi.5ap shows then that 
solutions of pi.Sp together with pl.7ap are also solutions of pi. 51) and viceversa. As a final remark, we recall that in 
O'Hanlon theory the standard conservation laws V'Sij — hold. An explicit proof of the vanishing of the covariant 
divergence of the energy- momentum tensor in modified theories of gravity can be found in (20| . 



B. The Cauchy problem for O'Hanlon gravity 

Taking into account the above dynamical equivalence, the Cauchy problem for f{R)- gravity can be defined as 
the Cauchy problem for the corresponding O'Hanlon theory. In this perspective, we discuss the Cauchy problem for 
the O'Hanlon gravity. Let us show the well-posedness of the Cauchy problem for system pi.5a|) and pi.6[) in vacuo. 
As we shall see, the same conclusions hold in presence of matter sources satisfying the standard conservation laws 
V*E,;j = 0. To this aim, we use generalized harmonic coordinates, given by the conditions 



F' - i7 ' = 



with 



P9' 



-V'lp 



(11.9) 
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The generalized harmonic gauge pi.9[) is a particular case of the one introduced in J27] to prove the well-posedness of 
the Cauchy problem for more general scalar-tensor theories of gravity As we shall see, the gauge pi.9p allows us to 
develop a second order analysis very similar to the one used in GR |23|. We rewrite eqs. (jll.SaP in the form 



Rij — 



where 



plays the role of effective energy-momentum tensor. The Ricci tensor can be expressed as [22 

1 



Ri 



with 



[g.,d,{F^ + HP)+g,,d.{F^ + H^)] 



= -^gP''d'pqg,,+A,{g,dg) 



(11.10) 



(11.11) 



(11.12) 



(11.13) 



where only first order derivatives appear in the functions Aij . Assuming ^ and taking the expression of _ff * into 
account, we obtain the following representation 



R^ 



+ -d^jV + B^j (.9, ip, dg, dip) 
ip ■' 



(11.14) 



where the functions Bij depend on the metric the scalar field p> and their first order derivatives. Analogously, using 
eq. pi.6p to replace all terms depending on the divergence g^'^Vp^ qtp, the right hand side of pi. 101) can be expressed 
as 



1 

^ L 



1. 



dlv + CtA9^^^dg,d(p) 



(11.15) 



Again, in the functions C'ij, only first order derivatives are involved. A direct comparison of eqs. pi. 141) with eqs. 
(III.15|) shows that, in the considered gauge, eqs. (III. 91) assume the form 



9^''9pggtj = Dij{g, ip, dg, dip) 



(11.16) 



Eqs. pi.l6p . together with eq. pi.l4p . form a quasi-dia gon al, quasi-linear second-order system of partial differential 
equations, for which well known theorems by Leray [l9l . |22| . [2a | hold. Given initial data on a space-like surface, the 
associated Cauchy problem is then well-posed in suitable Sobolev spaces [22||. Of course, the initial data have to 
satisfy the gauge conditions F^ = as well as the Hamiltonian and momentum constraints 



G' 



-r 



•Oi 



0,...,3 



(11.17) 



on the initial space- like surface. In connection with this, we notice that, from eq. pi.6p . we can derive the expression 
of the second partial derivative dgip and replace it in the right hand side of pi. 171) . so obtaining constraints involving 
no higher than first order partial derivatives with respect to the time variable . To conclude, we have to prove that 
the gauge conditions F^ = are preserved in a neighborhood of the initial space-like surface. To this end, we first 
verify that the divergence of the Einstein-like equations pl.Sa^ vanishes, namely 



V {ipG.j ~ T,j) = 

A straightforward calculation yields 

V' (^Q, - T,,) = (W) R.J - - ^) + ^^'G-,,- 



(11.18) 



(11.19) 
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By definition, the Einstein and Ricci tensors satisfy the identities V^Gij = and (VV) F^ij = (V'ViVj — VjV*Vi) (p. 

On the other hand, R — = is assured by field equations (|II.5bp . Therefore, identities (jll.lSP follow. If now gij 

dip 

and (f are the solutions of reduced Einstein-like equations (|II.16p and field equation (|II.6p . one has 

^G^' - T^' = - f (g'^dpF^ + g^PdpP; - g^^dpFP) (11.20) 

Identity (jll.lSP shows then that the functions satisfy necessarily the linear homogeneous system of wave equations 

gP'^dl^F; + E;W,FP = (11.21) 

where E^"^ are known functions on the space-time. Since the constraints (IIL17P amount to the condition doF^ (2^ 
on the initial space-like surface, a well known uniqueness theorem for differential systems like (IIL21|) (2^ assures that 
= in the region where solutions of pL16[) and (III.6I) exist. 
As mentioned above, the well-posedness of the Cauchy problem can be proved also in presence of coupling with 
standard matter sources, such as electromagnetic or Yang-Mills fields, (charged) perfect fiuid, (charged) dust, Klein- 
Gordon scalar fields. When this is the case, eqs. (IIL6I) and (|II.16p have to be coupled with the matter field equations. 
Applying the same arguments developed for GR [22j, .24-26], it is easily seen that, in the generalized harmonic gauge 
pi.9p . the matter field equations together with eqs. ()II.6[) and (jll.lGP form a Leray hyperbolic and causal differential 
system admitting a well-posed Cauchy problem. In addition to the well-known results by Bruhat's, the key point is 
that the field equations of matter field imply the standard conservation laws V'Sy- — 0. This fact allows to verify the 
validity of eqs. (jll.lSp in presence of matter too. 



III. METRIC- AFFINE /(ii)-GRAVITY 
A. Generalities 



In the metric-affine approach to /(i?)-gravity, the (gravitational) dynamical fields are a pseudo-Riemannian metric g 
and a linear connection T on the space-time manifold M . Specifically, in the so-called Palatini approach, the connection 
r is torsionless and it is not requested to be metric-compatible. In the approach with torsion, the dynamical connection 
r is forced to be metric but it has torsion different from zero. The field equations are derived from an action functional 
of the form 



-4(5, r)- I (vl5l/(^) + An) ds (III.l) 



where now R{g, T) = g^^Rij (with Rij := R'^^i^j) is the scalar curvature associated with the connection F. Throughout 
the paper we shall assume that the material Lagrangian is independent of the dynamical connection. In this case, the 
field equations are 



/'(«)«« - Tif {1)9.1 = S« , {III.2a) 



for /(i?)-gravity with torsion |27H29l |. and 



/'(i?)i?y - i/(i?)g., - , (111.3a) 



Vfe(/'(i?)g„)=0, (III.3b) 



for /(_R)-gravity in the Palatini approach |3ll - [33| . Considering the trace of Eqs. pil.2ap and pil.3a[) . we get the 
relation 



f'iR)R~2f{R) = i: 



(III.4) 
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which hnks the curvature scalar R with the trace of the stress-energy tensor S := g^^T,ij. From now on, we shall 
assume that the relation piL4[) is invertible that is E ^ const, (this implies, for example, f{R) ^ aR^ which is only 
compatible with E = 0) . Under these hypotheses, the curvature scalar R can be expressed as a suitable function of 
E, namely 

R = F{T.) . (III.5) 
If E = const., GR plus the cosmological constant is recovered [23|. Defining the scalar field 

^ /'(^^(E)) (III.6) 



we can put the Einstein-like field equations for both d la Palatini and with torsion theories in the same form |27l |33| , 

that is 

p 1 ~ _ 1 1 / 3 dip dip ^ 9(p dip 1^1. 

~ dip ^ ^^^^-^^ 

where we have introduced the effective potential 



V{ip) := - [ipF-\if')-\ip)) + ip'ifT'if)] 



(111.. 



for the scalar field if. In Eq. piI.Tp . Rij, R and V denote respectively the Ricci tensor, the scalar curvature 
and the covariant derivative associated with the Levi-Civita connection of the dynamical metric g^j. Therefore, if 
the dynamical connection F is not coupled with matter, both the theories (with torsion and Palatini-like) generate 
identical Einstein-like field equations. 

In addition to this, it has been shown that the Einstein-like equations (jlll.Tp (together with eas. (|III.6p ) are 
deducible from a scalar-tensor theory with Brans-Dicke parameter ljq — —3/2. To see this point, we recall that the 
action functional of a (purely metric) scalar-tensor theory is given by 



where ip is the scalar field, ip^ 



dip 



Wo 



ipup"- -U[ip)] + C. 



ds 



(III.9) 



and U{ip) is the potential of ip. The matter Lagrangian Cm{gij,'>P) is a function 



dx^ 

of metric and some matter fields ijj] ojq is the so called Brans-Dicke parameter. 

The variation of (jIII.9p with respect to the metric and the scalar field yields the field equations 



1 ~ _ 1 

Rij — 'l^RQij — 



and 



Wo 



1 



R 



ip'^ 



U' = 



(III.IO) 



(III.ll) 



where E, 



one obtains the equation 



and U' 



dU 

dip 



Taking the trace of eq. (jIII.10[) and using it to replace R in eq. pil.lll) . 



(2wo + 3) Vhip'' ^T. + ipU' -2U 



(III.12) 



3 2 

By a direct comparison, it is immediately seen that for wo •= and U(ip) := —V{ip) (where V(ip) is defined in eq. 

2 ip 

(jlll.Sp ) eqs. (jlll.lOp become formally identical to the Einstein-like equations (|III.7P for a metric-affine f{R) theory. 
Moreover, in such a circumstance, eq. piI.12p reduces to the algebraic equation 

6_ 



E + 2V'{^) - -V{^) 







(III.13) 
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relating the matter trace S to the scalar field ip, exactly as it happens for metric-affine /(i?)-gravity. In particular, it 
is a straightforward matter to verify that (under the condition /" ^ [27|) eq. (|III.13p expresses exactly the inverse 
relation of (jIII.6|) . namely 

Y. + 2V'{ip)~-V{^)=Q ^ J: = F-\{r)-\^)) (111.14) 

being F^^{X) = f'(X)X — 2f{X). When the dynamical connection does not couple with matter, metric-affine f{R)- 
gravity (with torsion or Palatini-like) are then dynamically equivalent to scalar-tensor theories with a Brans-Dicke 
3 

parameter uiq = — — . 

Moreover, we notice that field equations piL7l) can be simplified by passing from the Jordan to the Einstein frame. 
Indeed, performing the conformal transformation gij — (pgij, eqs. (|III.7p assume the simpler form (see for example 

MM) 

- 1%^. = Ie,^. _ ^V{cp)g,j (III.15) 

where Rij and R are respectively the Ricci tensor and the curvature scalar derived from the conformal metric gij . 

In connection with this fact, we recall that the conservation laws existing in the Jordan and in the Einstein frame, 
as well as their relations, are given by the following [s^, [Hi 

Proposition III.l Equations pil.7p . piI.SP and piI.13|) imply the standard conservation laws V-'Ey = 0. 

Proposition III. 2 The condition V-'Ey ^ is equivalent to the condition V^Tij = 0, where 
Tij :— —^ij ^V{(p)gij and V denotes the covariant derivative associated to the metric gij. 



B. The Cauchy problem in metric-affine formulation 

As mentioned above, if the trace of the stress-energy tensor E = const, both Palatini-like and with torsion theories 
reduce to GR with cosmological constant. Therefore, in such a circumstance, the Cauchy problem is well- formulated 
and well-posed [l^ [22[, l25j . For example, this is the case in vacuo and in presence of electromagnetic (or also Yang- 
Mills) fields {ii f{R)^aR^). 

On the contrary, if S ^ const., the situation is more complicated. In this case, the theory does not longer reduce 
to GR and so the well-formulation and well-posedness of the Cauchy problem is not automatically assured by the 
classical Bruhat results. 

In analogy with the purely metric case, matter sources could be again considered for Cauchy problem by using 

3 

the dynamical equivalence with scalar-tensor theories adopting a Brans-Dicke parameter ujq — — — . Un fortunately, 

here a difficulty occurs. The point is that the d'Alembertian gP'^WpWqLp disappears from eq. piI.12[) . It follows 
the impossibility of deriving the expression of the d'Alembertian g^'^'S/ q^p as a function of the dynamical variables 
and their derivatives up to the first-order. In other words, we cannot eliminate the second-order derivatives of the 
scalar field <p from the Einstein-like equations (jlll.lOp . This fact led some authors [33| to the conclusion that Palatini 
/(i?)-gravity has an ill-formulated Cauchy problem, so claiming a sort of no-go theorem for these theories. 

In the subsequent discussion, we shall prove that the conclusions achieved in |,37,] are not always true. Indeed, 
using the conformal transformation technique to pass from the Jordan to the Einstein frame, we shall give sufficient 
conditions ensuring the well-posedness of the Cauchy problem for metric-affine /(_R)-gravity coupled with a perfect 
fluid ^34. ,38i] or a Klein-Gordon scalar field [36]. The stated conditions will result as suitable requirements imposed on 
the function f{R), so representing a sort of selection rules for viable f{R) models. Moreover, we shall show that the 
set of viable f{R) models is not empty since the function f{R) = R + aR^ satisfies the above mentioned conditions. 

C. The Cauchy problem in presence of a perfect fluid 

Let us discuss the Cauchy problem for metric-affine /(i?)-gravity coupled with a perfect fluid. It is possible to show 
that, passing to the Einstein frame through a conformal transformation, the analysis of the initial value problem can 
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be carried out by following the same arguments developed in [23, [l^, [2^ for GR. Let us consider the metric gij of 
signature ( — h ++) in the Jordan frame and a perfect fluid with stress-energy tensor of the form 

Sii = (p + p) Ui Uj + p g,j (III. 16a) 

The corresponding matter field equations are 

VjS'^' = (111.16b) 

In eqs. (jIII.16[) . the scalars p and p denote respectively the matter-energy density and the pressure of the fluid, 
while Ui indicate the four velocity of the fluid satisfying the condition g'^WJJj = — 1. If we perform the conformal 
transformation gij = V'Sijj can express the field equations in the Einstein frame as 



- T;R9ri ^ (111.17a) 



and 



where 



VjT'J = (III. 17b) 

T., = i(p + p)C/,t/,+ (^-^) -g., (111.18) 

plays the role of the effective stress-energy tensor. Due to Proposition IIII.21 eqs. (|III.17bp are equivalent to eqs. 
(jlll.lGbl) . This point is crucial in our discussion, because it allows us to apply to the present case the results achieved 
in [22, 24, 25]. For the sake of simplicity, we shall suppose that the scalar field (p is positive, that is > 0. The 
opposite case ip < differs from the former only for some technical aspects and it will be briefly outlined after. Of 
course, using the above conformal transformation, it is implicitly assumed that 7^ at least in a neighbourhood of 
the initial space-like surface. 

Under these assumptions, the four velocity of the fluid in the Einstein frame can be expressed as Ui — ^/(pUi, so 
that the stress-energy tensor (jlll.lSP can be rewritten in terms of the four velocity Ui as 

T,, = ^{p+p) + - ^) -9., (111.19) 

After that, let us introduce the effective mass-energy density 



p- := 4 -t- ^ (III.20a) 



and the effective pressure 



P- ^ - ^ (III.20b) 



ip^ ip'-' 

In view of eqs. (jIII.20[) . the stress-energy tensor piI.lQp assumes the usual form 

T,j ^{p + p) tj.tjj + p-g,j (III.21) 

It is worth noticing that, starting from an equation of state of the form p = p{p) and assumin g that relation pil.20b|) 
is invertible (p = p{p)), by composition with eq. (jIII.20ap . we derive an effective equation of state p = p{p). In 
addition, we recall that the explicit expression of the scalar field (p as well as of the potential V{tp) are directly related 
with the particular form of the function /(i?). Then, the requirement of invertibility of the relation pil.20bp together 
with the condition ip> > (or, equivalently, ip> < 0) become criteria for the viability of the functions f{R). In other 
words, they provide us with selection rules for the admissible functions f{R) (see also (39|V 

At this point, the discussion of the Cauchy problem can proceed following step by step Bruhat's arguments (22l . 
[13, [25|. Bruhat's results are well known in literature, so we do not repeat the analysis here. We only recall that the 
Cauchy problem for the system of differential equations piI.17p . with stress-energy tensor given by eq. pil.2ip and 
equation of state p = p{p), is well-posed if the condition 
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is satisfied. It is worth noticing that, in order to check the requirement pil.22[) . it is not necessary to invert explicitly 
the relation pil.20bp . but simply we can verify the relation 

f = P^>l (111.23) 
dp dp /dp 

directly from expressions (jIII.20|) and equation of state p — p{p). Once again, condition (|III.23P depends on the 
expressions of (p and V{(p), therefore it is strictly related to the form of the function f{R). Thus, condition pil.23[) 
represents a further criterion for the admissibility of /(_R)-models. 

For the sake of completeness, we outline the case cp < 0. If we still suppose that the signature of the metric in the 

Jordan frame is ( — h ++), the signature of the conformal metric is (H ) and the four velocity of the fluid in the 

Einstein frame is then Ui = —ipUi- The effective stress-energy tensor is given now by 

T^, = -^(p + p) US J + - 9^3 = (P + P) - Pmj (III.24) 

where we have introduced the quantities 



and 



p^ (p3 



(III.25a) 



(III.25b) 



representing, as above, the effective mass-energy and the effective pressure. At this point, everything proceeds again 
as in (2^ . [23 , [2^, except for a technical aspect: now the quantity r := p + p = ^ is negative (if, as usual, p 

and p are assumed positive). In this regard, the reader can easily verify that, with the choice log(— /^^r) instead of 
log(/^'^r) as in (2^. l24l [25|. the Bruhat's arguments apply equally well. 

Let us now give a simple and illustrative example of the above results. Let us take into account the /(i?) = R + aR^ 
model coupled with dust (for, example, let us consider a cosmological model). The matter stress-energy tensor in the 
Jordan frame is — pUiUj. The trace of the Einstein-like equations pil.2a[) yields the relation 

(1 + 2aR)R -2R- 2aR^ = -p i — > R^ p (III.26) 
The scalar field pil.6p is given by 

ifip) = f'{R{p)) = l + 2ap (III.27) 

For small values of the density p << 1 (for example, the present cosmological density) and values of \a\ not comparable 
with 1/p, we can reasonably suppose p > 0, independently of the sign of the parameter a. Let us calculate now the 
potential piI.SP 

vip) = i [pF-\ifr'm + ^'{fTH^)] (ni.28) 

Since {f')^^{p) = p, from eq. pil.27p we have 

i^2(/')-^(^) = i(l + 2ap)V, (111.29) 
also, considering the relation F-^{Y) = f'{Y)Y - 2f{Y), it is 



and 



^F-\{,frH^)) - \f-\p) = -p. (IIL30) 



l^F-\ifr\^))^-i^±f^, (III.31) 
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We can conclude that 



VMp)) - '^P^:^ (111.32) 
In the Einstein frame, the stess-energy tensor is expressed as (comparing with eq. pil.191) ) 

T,, = A^U.U, - ^-g., (III.33) 
The latter can be seen as the stress-energy tensor of a perfect fluid with density and pressure given respectively by 



and 

V{^) ap" 



(p3 2(1 + 2apY 

Under the above assumptions, it is easily seen that function pil.34bp is invertible. For p > 0, we have 

dp Aap 



(III.34b) 



dp 4(1 + 2apf 

Moreover, we have 

dp 4 — Aap 



dp 4(1 + 2apY 

so that 

dp dp I dp —1 + ap 



^ (III.35) 



(III.36) 



dp dp /dp ap 



> 1 ^ a<0 (III.37) 



The conclusion is that the model f{R) = R + aR^, with a < 0, has a well-posed Cauchy problem when coupled with 
dust. This is in agreement with results in (40j . 



D. The Cauchy problem in presence of a scalar field 

Let us consider now the Cauchy problem for metric-afhne /(i?)-gravity coupled with a Klein-Gordon scalar field. 
As above, we shall state sufficient conditions ensuring the well-posedness of the problem. Let -0 be a Klein-Gordon 

scalar field with self-interacting potential U{4>) = -rri^ip'^. The associated stress-energy tensor is expressed as 



The corresponding Klein-Gordon equation is given by 

V,^5'^=mV (nL39) 
where V denotes the Levi-Civita covariant derivative induced by the metric g^ . The trace of tensor (|III.38p is 

E:=E,-,-.-|^^,--2^V. (nL40) 

We notice that trace pil.401) depends explicitly on the metric tensor gij. In view of this fact, we cannot apply the 
conformal transformation directly to the field equations (IIII.7P , being the scalar field ip defined by pil.6p . Indeed, if 
we proceed in this way, both the metric gij and g^j would appear in the conformally transformed equations piI.lSp . 
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3 

The idea to overcome this difficuhy is using the aheady mentioned dynamical equivalence with ujq = —— Brans-Dicke 

3 

gravity. In other words, we discuss the Cauchy problem for a loq = Brans-Dicke theory coupled with the given 

Klein-Gordon field ip. The field equations of such a theory are the Einstein-like equations pil.7|) . the equation piI.13p 
and the Klein-Gordon equation pil.39|) . We perform the conformal transformation gij — (pgij. The scalar field ip is 
not longer defined by piL6p , but it is now a dynamical variable related to the trace S through eq. (IIIL13P . This is an 
important difference with respect to the perfect fluid case discussed above. Anyway, after conformal transformation, 
the Einstein-like equations pil.7|) assume the simpler form pIL15|l . Moreover, using the relation 



f - - f - -'^ + —^s'^ - — ^af\ + —^s^ an 41) 



between the Levi-Civita connections associated respectively to gij (Tij^) and gij (Fj^'*), it is easily seen that the 
Klein-Gordon equation expressed in terms of the conformal metric gij becomes 

where Vj denotes the covariant derivative associated to the conformal metric gij. Also the trace E can be expressed 
as function of gij , that is 

dip 

Now, the relation pil.isp links the scalar field ip to the Klein-Gordon field ip, its partial derivatives — — and the 

. ox^ 

conformal metric gij. As already said, in eqs. piI.15p the quantity 

T.J := ^S., - ^Viv)9^J (111.44) 

plays the role of effective stress-energy tensor. The Klein-Gordon equation pil.39[) implies the conservation laws 
V-'Eij = and thus also the identities V^Tij = (due to Proposition IIII.2I) . The latter is a crucial fact which allows 
to use harmonic coordinates, once again following the same arguments developed in (l9l. [^S. [25j . 
More in detail, we rewrite the Einstein-like equations pil.isp in the equivalent form 



% = T^J - T^Tg.j (III.45) 



1, 

Adopting harmonic coordinates, i.e. local coordinates satisfying the condition 

VpW = -5P«f;^ = 0, (III.46) 
we can express equations pil.45[) as (see, for example, (l9l. [25|l 

.V, ^^9^J -f^^(^-g^d-g,i^,d^) (III.47) 



dxPdxi 



where indicate suitable functions depending only on the metric the scalar field ip and their first order derivatives. 

Moreover, we assume that equation piI.13|) is solvable with respect to the variable (p. In other words, we suppose 
to be able to derive from eq. piI.13p a function of the form 

expressing the scalar field ip in terms of the metric g, the Klein-Gordon field ip and its first order derivatives. In 
particular, from eq. (1111.43^ . it is easily seen that the dependence of (p on the derivatives of ip is necessarily of the 
indicated form. The requirement about equation pII.13P depends on the explicit form of the potential V{(p). The 
latter is determined by the function f{R) through the formula pil.81) . Therefore, the above assumption becomes a 
rule to select viable /(i?)-models. In addition to this, from eq. (IIII.48|) . we derive the identity 



^/ dtp dip 
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Inserting eq. pil.49p in eq. (jIII.42[) and taking eqs. (IIII.46[) into account, we obtain the final form of the Klein-Gordon 
equation given by 

In eqs. pil.49[) and pil.50p . fi and / indicate suitable functions of gij, ip and their first order derivatives only. 

Eqs. pil.47p and pil.50p describe a second order quasi-diagonal system of partial differential equations for the 
unknowns gij and ip. The matrix of the principal parts of such a system is diagonal and its elements are the 
differential operators 

gP'^ ^ (111.51a) 

and 

' r - - . . ^r^r\ ^ (ni.5ib) 



\ ax^ ox^ ^ 



The operator pil.Slap is noting else but the wave operator associated with the metric gij , while the operator pil.51b[) 
is very similar to the sound wave operator arising from the analysis of the Cauchy problem for GR coupled with 
an irrotational perfect fluid (22l . [25| . To discuss the Cauchy problem for system pil.47[) and pil.50[) . we can then 
follow the same arguments developed in [2^ [2^. In particular, if the quadratic form associated with pil.51b[) is 
of Lorentzian signature and if the characteristic cone of the ciperator pil.51bp is exterior to the metric cone, the 
system pil.47|) and (IIII.SOP is causal and Leray hyperbolic (23i. l2a| . In such a circumstance, the corresponding Cauchy 
problem is well-posed in suitable Sobolev spaces. Still following [2^. l25j. if the signature of gij is (-1 ), the 

required conditions are satisfied when the vector — — -g^-' is timelike and the inequality 



dxi' 

dip 



> (III.52) 



holds. On the contrary, if the signature of gij is ( — h ++), the inequality pil.52p has to be inverted. As already 
mentioned, the explicit expression of the function pll.48|) depends on that of the potential pll.8|) which is determined 
by the function f{R). Therefore, the requirement (1111.52^ (or, equivalently, its opposite) can be a criterion to single 
out viable /(_R)-models. 

As an example, let us consider again the model f{R) = R + aR^. From the identities F^^{X) = f'{X)X — 2f{X) = 
—X and {f')^^{(p) = ^- — as well as the definition pll.81) . we get the expression of the effective potential 

l/(^) = -^(^-l)V (111.53) 
Eq. pll.53|) together with eqs. pil.isp and pil.43p yield 



(1 +2mV) 

ip = (III.54) 



f 2 ay- dip 

\ 2a dx^ dx^ 



gSt 



which expresses the scalar field ip as function of the metric gij, the Klein-Gordon field ^ and its first order derivatives. 
Directly from pil.54p it follows 

- (III.55) 



\dx^ dx^ J \2a dx^ dx^ 



In the signature (^ ) for the metric gij, it is immediately seen that the requirement pil.52p is automatically 

di/j 

satisfied if a < and if 9^'^'q~^ ^ time-like vector field. If the signature is ( — h ++), the condition becomes a > 0. 
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IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have discussed the initial value problem for /(i?)-gravity in metric and metric-affine approaches. 
Both the well-formulation and the well-posedness hold in metric formulation. This result is easily achieved as soon as 
the theory is recast in term of the O'Hanlon theory. In this case, by choosing a suitable gauge, the Bruhat arguments 
for GR can be easily reproduced for any form of matter acting as source (e.g. perfect fluid, electromagnetic, Klein- 
Gordon, Yang-Mills fields). In vacuum case, /(i?)-gravity is equivalent to GR plus cosmological constant. 

On the other hand, it can be shown that the initial value problem is, in general, well-formulated also for metric-afhne 
/(i?)-gravity, at least in presence of the standard matter sources mentioned above (35l. [39|. This means that there 
are no objections to the viability of metric-afhne /(i?)-gravity, based on the supposed ill- formulation of the Cauchy 
problem [? ]. However, also the well-posedness is necessary in order to achieve the complete control of dynamics but, 
in this case, the role of source fields and the specific form of the function f{R) have to be carefully discussed. 

At present, no general results exist about the well-posedness of the Cauchy problem for metric-afhne /(i?)-gravity. 
Therefore, to discuss the well-posedness, we are forced to proceed through model by model analysis. For example, 
anytime the trace of the stress-energy tensor is constant, the theory reduce to GR plus cosmological constant and the 
well-posedness is thus recovered, independently of the explicit form of the function f{R). In case of coupling with a 
perfect-fluid or a Klein-Gordon scalar field, sufhcient conditions can be stated in order to assure the well-posedness. It 
has been proved that there exist some functions f{R) that actually satisfy such conditions which then become a sort 
of selection rules to determine the viable forms of /(_R)-models. This feature is, in some sense, a sort of consistency 
check based on a first principle: the initial value problem selects viable models. From a physical viewpoint, this fact 
is crucial to discriminate among competing theories that have to be compared with data. 
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